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Abstract 

We examine the dynamics of biological nanomotors within a simple model of a rotor having three 
ion-binding sites. It is shown that in the presence of an external dc electric field in the plane of 
the rotor, the loading of the ion from the positive side of a membrane (rotor charging) provides 
a torque leading to the motor rotation. We derive equations for the proton populations of the 
sites and solve these equations numerically jointly with the Langevin-type equation for the rotor 
angle. Using parameters for biological systems, we demonstrate that the sequential loading and 
unloading of the sites lead to the unidirectional rotation of the motor. The previously unexplained 
phenomenon of fast direction-switching in the rotation of a bacterial flagellar motor can also be 
understood within our model. 



PACS numbers: 87. 16. Ac, 85.85.+j 



Biological rotary motors provide remarkable 



jEJxamples of how electrochemical energy can 
2] . Two of the most important representa- 



be efficiently converted into mechanical motion 
fives of this family, the F motor of ATP (adenosine triphosphate) synthase (FxFo-ATPase) 
and the bacterial flagellar motor (BFM), are powered by H + or Na + ions, which flow down the 
electrochemical gradient across the mitochondrial or cell membranes, thereby generating a 
torque 3|, |4| ■ Hereafter we concentrate on proton-driven motors, but the same considerations 
are valid for sodium-driven motors as well. The gradient of the electrochemical potential is 
maintained by the metabolic mechanism translocating protons from the negative side of the 
membrane to its positive side 

Both of these rotary motors have similar components: (i) a stator, tightly attached to 
the membrane, and (ii) a ring-shaped rotor, which can freely rotate around its axis. The 
rotor part of the Fo motor is mechanically coupled to the Fi domain of ATP synthase, 
whereas the rotor of the BFM is linked to the propeller. It is assumed j3] that the rotor 
has several (10 to 14) proton-binding sites. The permanent generation of the torque can be 
derived from the electrostatic interaction between stator charges and charges of the rotor 
sites which bias the thermal diffusion of the rotor in a specific direction jsl, 6|. This Brownian 
ratchet mechanism can generate a torque of about 40 pN-nm corresponding to a realistic 
rate of ATP synthesis [3]. However, both the power stroke and Brownian ratchet models, 
exploiting "constructive" features of Brownian motion, do not succeed in explaining the 
amazing performance of the bacterial flagellar motor, generating a torque of 2700 - 4600 
pN-nm with an efficiency of about 90% 4|. The ability of the BFM to rapidly switch the 
direction of the rotation remains unexplained as well, because transport in ratchets (see, 
e.g., the reviews in 7(]) is usually fully controlled by a fixed asymmetry of the potential 
energy. These facts point to the possibility that the BFM can use the energy stored in the 
proton electrochemical gradient directly, without a mediation of the Brownian motion (J. 

Here we explore a simple model mimicking important features of real biomolecular ro- 
tary motors and allowing a quantitative treatment based on methods of condensed matter 
physics 8]. These approaches have been previously applied to nanoelectromechanical sys- 
tems (NEMS) with their mechanical motion affecting the electrical properties of electronic 
devices |9|. Similar processes take place in nanoscale biological objects, where electrical 
and mechanical degrees of freedom are also strongly coupled, making them living counter- 
parts to artificial NEMS. Note that only nano-oscillators have been extensively studied by 



theorists, although a single- molecule rotor and a nanoelectromechanical rotational actuator 
lp| have been demonstrated experimentally. To the best of our knowledge, no theoretical 
investigations of rotary NEMS have been reported yet. 

Model. The system under consideration here consists of three equally-spaced proton- 
binding sites A, B, C attached to the rotating ring (rotor) in the presence of a constant y- 
directed electric force F (see the inset in Fig.l). The rotor sites can be coupled to two proton 
sources Si and £2 as well as to the proton drain D, connected to the proton reservoir with 
a low electrochemical potential \ir> (the negative side of the membrane). The source leads 
Si and 5*2 can be connected (or disconnected) at will to the proton reservoir with a higher 
electrochemical potential /is (the positive or P-side of the membrane). We show below that 
the activation of the lead Si results in a clockwise motion of the rotor, whereas connecting 
the SVlead to the P-side of the membrane (and disconnecting the Si-lead) generates a 
counterclockwise rotation. At each instant of time only one source lead is coupled to the 
P-side proton reservoir. 

The Hamiltonian of the system has the form: 



where <ft is the angle of rotation, p = —ih~(d/d<f)) is the operator of angular momentum of 
the rotator with radius r and effective mass M. We take into account here the effects of a 
constant y-directed external electric field, with a potential energy profile U(<f)) = — Uocosrf), 
on the protons localized in the three sites a = A, B, C with positions, characterized by the 
angles (ft a, (ft b, fie, respectively. The operators a+,a a describe the creation and annihilation 
of a proton on the site (dot) o with a population n a = a^a a , whereas the operators c£ a , c^ a 
are related to the /c-state of the proton in the source and drain reservoirs (leads) with energy 
Eka {a = Si, S2, D). The Coulomb repulsion between protons is given by the potentials U aa i. 
The tunneling coupling between dots and leads is given by the Hamiltonian 





(1) 




where the tunneling amplitudes T^a are multiplied by the factor 




which reflects an exponential dependence of the tunneling rate on the distance between 
the a— dot and the a— lead with a characteristic spatial scale A. To take into account the 
influence of the classical dissipative environment {Q} with the Hamiltonian .Heath on the 
rotational degrees of freedom, we include the term Hq = —r (j)Q + H Bath in Eq.(JTJ. Here 



rotor to the bath is characterized by the drag coefficient (. The unperturbed bath variables 
have Gaussian statistics with zero average, (Q^) = 0, and the correlation function, 
(Q( \t)Q(°> (t 1 )) = 2C ! T5{t — t'), where T is the temperature of the environment (k B = 1). 
The Brownian motion of the nanorotator is governed by the Langevin equation 



where 7 is the damping rate of the rotator, 7 = (/M, and the Gaussian fluctuation source 
f = Q(°) /(Afr ) is characterized by the correlation function: (f (t)f (*')) = ^(T/Mr^)6{t-f). 
The influence of the proton tunneling events, ~ —dH tun /d<ft, on the mechanical motion is 
assumed to be negligibly small here. 

To describe the process of loading and unloading of proton-binding sites A, B, and C, 
we introduce the proton vacuum state |1) = |Vac), jointly with seven additional states, 



|2> = a+|l>, |3> = a+|l>, |4> = a>+|l>, |5> = a+|l), |6) = a+a+|l>, |7> = a+a+|l>, |8> = 



a\a^a^\l). Each of the proton operators can be expressed in terms of operators p^ = 
\li){v\ (/i, v = 1,..,8). In particular, the operator a a has the form: a a = Ylav a ^ v P^i 
with the following non-zero matrix elements: a^ ; i2 = ^a-^a = QU;56 = a A ; 78 = 1;°B;13 = 
-a B . 2 A = a B . 57 = -a B . fiS = l;ac ; is = -acw = -ac ; 37 = ^c-as = 1- The populations 
of the dots, n a = a+a CT , are expressed in terms of the diagonal operators p M = as: 
riA = P2 + Pi + Pe + Ps, n B = p 3 + p 4 + pi + p 8 , n c = pn + Pe + Pi + Ps- Thus, for the protons 
localized on the rotor sites A, B,C we obtain the Hamiltonian Habc = Eu=i e M | /x) {/x | , 
with an energy spectrum depending on the local value of the rotor angle <fi: e% = 0, €2 = 
Ea-UqCOS^ + ^a), e 3 = E B -U cos((j) + (j) B ), e A = e 2 + e 3 , e 5 = E c - U cos(0 + <j> c ), e 6 = 
e 2 + e 5 , e 7 = e 3 + e 5 , e 8 = e 2 + e 3 + e 5 . We assume here that the characteristic time of proton 
tunneling to and out of the proton-binding sites, 7 _1 , is much shorter than the the time scale 
of the rotary angle, and that the noise produced by the proton tunneling between the 

sites and the source and drain contacts has much less effect on the mechanical motion of the 



Q is the bath variable, which can be represented as a sum of the fluctuating part, Q^°\ and 
the bath response on the action of the rotator, Q = — (ro<j), where the coupling of the 





rotor than the noise £ generated by the bath {Q}. Accordingly, we can average the stochastic 
Eq.(J2]) over fluctuations of the proton reservoirs without averaging over the fluctuations 
of the mechanical heat bath. The partially averaged proton population n a involved in 
Eq. (j2J) depends on the local fluctuating value of the rotational angle <fi. To determine these 
populations, we derive the following master equation for the proton distribution p M , averaged 
over reservoirs fluctuations, 

/V + 7 m Pm = with ( 3 ) 

V 

oca 

7 M = Jlulvm and T a (<f>) = T a \w a(7 ((j))\ 2 . Also: = e M -e^ {h = l),T a = 27r]T fc \T ka \ 2 5{u- 
uj^y) [8|. The protons in the reservoirs are characterized by the Fermi distributions, f a (uj) = 
[exp( ""^" ) + l] , with temperature T and electrochemical potentials = V/2,fi D = 
— V/2, where V is the proton voltage build-up. We include the absolute value of the proton 
charge, |e|, into the definition of the voltage V and measure the voltage in units of energy, 
meV. Notice that, despite of the averaging over proton reservoirs, the master equation, 
Eq.(J3]), contains a stochastic component, which is determined by the fluctuations of the 
rotor angle <j>(t), taken at the same time t. Eq.Q is valid for weak coupling between dots 
and leads, and for the case when the angular coordinate <p(t) does not change significantly 
on the characteristic time scale, tt = 7r/T, of the Fermi distribution. 

We consider here the rotation of the system with the parameters roughly corresponding 



111 ], which has mass M = 5000 kDa 



to the rotor part of the bacterial flagellar motor 
8.3 xlO -21 kg and a radius r = 15 nm. Three torque-generating sites, A, B, and C, are 
attached at the points 4>a = 0, <Pb — 2vr/3, and <pc — — 27r/3, respectively. The locations of 
the two possible source contacts Si, S2 and the drain D are defined as <j>s 1 = —2tt/3, <ps 2 — 
27r/3, 4>d = 0. For the Coulomb interaction between sites located a distance r \/3 apart, 
in a medium with a dielectric constant e r ~ 4, we obtain: U a b — Ub c — U ac ~ 15 meV. We 
choose an intermediate value of the drag coefficient ( = 30 pN s/m, which is related to the 
motion of a sphere with radius r = 15 nm in an environment with a viscosity of 0.1 mPa-s. 
The rotational damping rate is quite high, 7 = 3.6 ns -1 , which means that the free rotations 
of the motor come to an end after a time interval of the order 7 _1 ~ 300 ps. 

Results. In Fig. la we present a schematic diagram of the rotor ring, jointly with the time 



dependence of the number of full rotations 0(t)/27T, obtained from the numerical solution of 
the Langevin equation, Eq.(T5]), and the master equations, Eqs.®, at the source-drain voltage 
V = Hs—^d = 500 meV, the external potential U = 200 meV, the tunneling couplings to the 
leads = T^j = 10 9 s _1 , and at temperature T = 300 K. The source contact Si is activated, 
and, accordingly, the rotor moves in the clockwise (positive) direction and performs more 
than four full rotations in ~10 /is. We start here with a slightly shifted initial position of 
the rotor, when 0(0) = — n/20, 0(0) = 0, and the sites A and C are displaced from the 
drain D and source S\ contacts, respectively. The pronounced influence of the environment 
on the rotational degrees of freedom is reflected in the noisy time dependence (see Fig. lb) 
of the speed of rotations f2(t)/27r, where Q(t) = <fi(t). Figure lc illustrates the synchronous 
dynamics of loading and unloading the proton binding sites A, B, and C having populations 
riA, us-, and nc, when they pass through the source and the drain leads. The site C is 
populated first because of its close proximity to the source lead S\. The external electric 
field pushes the C-proton, and, correspondingly, the whole rotor unit, to turn through the 
angle 2tt/3 to the position of the minimum of the potential U = — Uocos(<f) + 0c*). At this 
position of the rotor, the site C approaches the drain contact D and unloads the proton. At 
the same rotor position, the site B starts to be populated since this site is in the loading range 
of the source lead S\. This leads to a subsequent 120°-turn of the rotor. The process repeats 
over and over, resulting in a continuous unidirectional rotation of the rotary ring. The 
torque, exerted by the nanomotor p, can be defined as: T = — J2 a (j 1 ^ ^^T^ / > where 
£7(0) = —Uq cos0 is the potential of the -(/-directed external electric field F = Fy, U = Fr . 
As follows from Fig. Id the torque oscillates between 10 and 30 pN-nm. These numbers are 
comparable with the torque produced by the Fo-motor of ATP synthase, but much less than 
the BFM torque. The torque of the motor can be increased by increasing the number of 
torque-generating sites. 

In Fig. 2 we present the dependence of the average speed of rotations, (Q)/2n = (<j))/2ir 
(left axis), and the average particle current, / = —Ii = —(d/dt)^2f.( c iL c kL) (right axis), 
on the difference of proton electrochemical potentials of the positive (/is) and the negative 
(Hd) side of the membrane (proton voltage build-up, V). At the fixed voltage V the average 
speed of rotation is linearly proportional to the amplitude of the external electric potential 
Uq, up to the threshold value U 0traaiX = V/2. 

To demonstrate that the unidirectional rotation of the motor has an origin which is 



different from the mechanisms utilizing Brownian motion of the rotary ring, we drop the 
fluctuation source £ in Eq. (T5]) and calculate the time dependence of the angle in the absence 
of noise (it corresponds to the mean-field approximation for rotational degrees of freedom). 
To show the reversibility of rotations, in Fig. 3 we switch also the direction of rotation from 
clockwise to counterclockwise by connecting the lead S2 to the proton reservoir with a higher 
electrochemical potential (and disconnecting the lead Si). The initial condition for the rotor 
angle remains the same, 0(0) = — 7r/20, 0(0) = 0. This shift of the rotor position hampers 
the initial loading of the site B from the source S% (see Fig. 3c), postponing the beginning of 
the full-scale revolutions. Nevertheless, without noise, the rotor makes more than three full 
rotations in 10 /is, despite the unfavorable initial conditions and despite the strong enough 
damping rate 7 = 3.6 ns -1 . Note that the time dependence of the instantaneous speed 
of rotations Q/2tt (see Fig. 3b) shows regular oscillations, compared with the noisy time 
dependence of the same variable depicted in Fig. lb. It is evident from Fig.3b that the speed 
of rotations oscillates synchronously with the torque generated by the motor (see Fig. 3d). 

Turning back to the situation with noise, we now estimate the efficiency of the nanomotor 
proposed above. For the parameters Uq = 200 meV, and V = 500 meV, the average torque 
(T) is about 20 pN-nm, while the average speed of rotations (Q)/2tt is near 420 kHz (see 
Fig. Id and Fig. 2). This means that for the output power of the motor Q] we obtain: 
(T) • (O) ~ 3.3 • 10 5 eV/s. To calculate the input power provided by the proton reservoirs, 
we assume that the proton flux (proton current) is about I = 1.2 • 10 6 particles per second 
(see Fig. 2), and the voltage drop between the reservoirs is V — 0.5 eV. Thus, the input 
power can be evaluated as I • V = 6 • 10 5 eV/s. Therefore, for the efficiency of the motor, 
Eff = (T) • • V), we obtain the estimate: Eff ~ 55%. Increasing the number of 

protonable sites and optimizing the system's design and parameters would improve the 
performance of our simple model, approaching the performance of real rotary biomotors, 
having an efficiency near 90%. 

In summary, we have proposed a simple physical model of ion-driven rotary nanomotors 
which transmit the energy of an external constant electric field into mechanical motion in the 
presence of an electrochemical gradient. This model describes the main properties of both the 
Fo-motor of ATPase in mitochondria membranes and the bacterial flagellar motor without 
resorting to Brownian motion. The proposed mechanism of torque generation is based on 
the subsessive loading and unloading of ion-binding sites on the rotor in the presence of a dc 



electric field. Our model can explain the previously mysterious phenomenon of rotational 
direction switching which occurs in bacterial flagellar motors. Solving jointly the Langevin- 
type equation and master equations for the populations of the ion-binding sites, we have 
determined the time evolution of the system for various values of voltages and amplitudes 
of the dc external electric field. 
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FIG. 1: (Color online) (a) Time dependence of the number of rotations (j){t)/2ii at V = 500 
meV, Uq = 200 meV, and T = 300 K, with a schematic diagram of the system in the inset; (b) 
Stochastic dynamics of the instantaneous speed of rotation, 0,(t)/2n = c/)(t)/2ir; (c) Populations of 
the proton-binding sites versus time; (d) Time dependence of the rotor torque in the presence of a 
heat bath. Notice the periodicity in (c) and (d). 
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FIG. 2: (Color online) Dependence of the average speed of rotation, {Q)/2tt (solid blue curve, left 
axis), and the average proton current, / (dashed black curve, right axis), on the proton voltage 
build-up, V = ld S - fiD, at U = 200 meV and T = 300 K. 




FIG. 3: (Color online) Time evolution of the system in the absence of external noise at V = 500 
meV and Uq = 200 meV, but now switching the direction of rotations: (a) number of full rotations, 
(j>(t)/2ir, as a function of time; (b) speed of rotations, Q(t)/2ir; (c) populations nA(t),riB(t), and 
nc(t), describing the loading and unloading of the sites; (d) torque exerted by the rotor versus 
time. 



